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We investigate intermittent difltusion using cycle expan- 
sions, and show that a truncation based on cycle stability 
achieves reasonable convergence. 

PACS: 5.45.+b 



I. INTRODUCTION 

Classical dynamical systems range from purely inte- 
grable to purely hyperbolic. For purely integrable sys- 
tems we have a variety of classical methods, such as sep- 
aration of the Hamilton- Jacobi equation [h|. For almost 
integrable systems we have KAM theory M. For purely 
hyperbolic systems it is possible to obtain much infor- 
mation about the system by grouping contributions com- 
puted on unstable periodic orbits |^ into terms in cycle 
expansions |^J^. They yield the classical escape rate of 
open billiard systems to a high degree of accuracy [||, 
and the semiclassical energy levels of systems such as he- 
lium [0 using a surprisingly small number of unstable 
periodc orbits. 

However, the formalism does not work well for generic 
dynamical flows for which the hyperbolic regions coexist 
with attractors, intermittent regions and elliptic regions. 
For intermittent systems the cycle expansions ordered 
by the topological cycle length converge poorly if very 
long almost stable cycles dominate the dynamics. The 
original Q as well as more recent applications [|| of cycle 
expansions to intermittent systems use detailed analytic 
information about the intermittent regions in order to 
explicitly sum infinite sequences of such cycles. 

Our philosophy here is that it should be possible to ob- 
tain reliable dynamical averages without a complete un- 
derstanding of the detailed structure of the phase space, 
as long as we are restricted to a given connected region 
in which the dynamics is ergodic. Recent work |^ on the 
Lorentz gas suggests that reordering the cycle expansions 
by stability pO[ may improve convergence in such situa- 
tions. Here we test this proposal by calculating diffusion 
in a one-dimensional intermittent map and demonstrate 
that the stability ordering yields better convergence than 
the ordering by the topological cycle length. 

There are several arguments in favour of using stabil- 
ity rather than the topological or (in case of continuous 



flows) real time length as the truncation criterion: 

1. Longer but less unstable cycles can give larger 
contributions to a cycle expansion than short but 
highly unstable cycles. In such situation trunca- 
tion by length may require an exponentially large 
number of very unstable cycles before a significant 
longer cycle is first included in the expansion. 

2. Stability truncation requires only that all cycles up 
to given stability cutoff be determined, without re- 
quiring detailed understanding of the topology of 
the flow and symbolic dynamics. It is thus much 
easier to implement for a generic dynamical sys- 
tem than the curvature expansions H| which rely 
on finite subshift approximations to a given flow. 

3. The stability ordering preserves approximately any 
shadowing that is present. That is, a long cycle 
which is shadowed by several shorter ones will have 
a stability eigenvalue which is approximately the 
product of the shorter cycle eigenvalues, and will 
be most likely be included at the same stability 
cutoff. 

4. Cycles can be detected numerically by searching a 
long trajectory for near recurrences [ pTp^ . The 
method preferentially finds the least unstable cy- 
cles, regardless of their topological length. Another 
practical advantage of the method (in contrast to 
the Newton method searches) is that it only finds 
cycles in a given connected ergodic component of 
phase space, even if isolated cycles or other ergodic 
regions exist elsewhere in the phase space. 

In what follows we illustrate the first three points by 
investigating the convergence of stability cutoff approach 
for a simple system. We begin by describing diffusion 
on a lattice of one-dimensional maps, how to calculate 
the diffusion coefficient using cycle expansions, and then 
perform the calculations numerically. Finally we discuss 
the scope of such approaches and possible improvements. 

II. DIFFUSION IN 1-D MAPS 

As a model on which to test the above ideas we shall 
use a well understood one parameter family of diffusive 
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one-dimensional maps. For such maps the symbohc dy- 
namics is a complete binary shift, all cycles can be ex- 
haustively enumerated, and the limitations of the length 
truncated cycle expansions are solely due to the lack of 
hyperbolicity, and not to inadequate understanding of 
the symbolic dynamics. 

Many of the features of intermittent systems can be 
captured by one-dimensional intermittent maps intro- 
duced in ref. Q to study turbulence. Piecewise lin- 
ear approximations [p^ can make statistical mechanics 
aspects of such intermittent dynamics, including phase 
transitions and Levy flights, analytically tractable. 
Intermittent maps can lead to anomalous deterministic 
diffusion, with the mean square displacement either sub- 
linear or superlinear in the time 16|. 

In the interval x G [—1/2, 1/2), which we call the ele- 
mentary cell, our model map takes the form 



X = X 



[x + 1/2] , 



/(x) =^(l + 2|2x|") 



(1) 



where a > —1. For any value of a, this maps the interval 
monotonically to [—3/2,3/2). Outside the elementary 
cell, the map is defined to have a discrete translational 
symmetry. 



f{x + n) = f{x) + n 



ne Z 



A typical initial x in the elementary cell diffuses, wan- 
dering over the real line. The map is parity symmetric, 
f{—x) = — /(i) , so the average value of — Xn is 
zero, and there is no mean drift. 




FIG. 1. The map (|l|) for a — 1 restricted to the elemen- 
tary cell, together with the 0®+ cycle. In the full space this 
cycle correspond to ballistic motion to the right, with each 
7th iteration entering the next cell. 

We now restrict the dynamics to the elementary cell, 
that is, we define 



where [z\ is the greatest integer less than or equal to z, so 
that X is restricted to the range [—1/2, 1/2). The reduced 
map (see Fig. ||) is 



f{x) - f{x) - [f{x) + 1/2] 



(2) 



Xj with sta- 



A cycle p = {xi,X2, ■ ■ ■ ,Xn}, P{x-j) = 
bility Ap = 11^=1 /'(^^i) corresponds to a trajectory 
which returns to an equivalent point in the full dynamics, 



fip. Thus the cycles fit into two categories, 



those which are periodic in the full dynamics, fip 



0, 



and those which are not. The diffusive properties of the 
map are fully specified by the reduced map /(x), together 
with the lattice translation n. The diffusion constant 



D ^ lim 



(3) 



is computed as an average over initial conditions in the 
elementary cell. 

The reduced map (^) has three branches, correspond- 
ing to either moving to the left, staying in the elemen- 
tary cell, or moving to the right, h = { — 1,0,1}. Hence 
the natural symbolic dynamics is a 3-letter alphabet 
{—,0,-1-}. For a given symbol string the total transla- 
tion n is just the sum of the individual symbols in the 
cycle symbol string. The three branches form a Markov 
partition, because each is mapped onto the whole inter- 
val, and the symbolic dynamics is thus unrestricted in 
the three symbols, with all finite strings corresponding 
to cycles. 

The point a; = is a fixed point (cycle of length 1) with 
symbol sequence 0. For a < 0, Aq = /'(O) = oo, this 
fixed point is infinitely unstable, and its contribution to 
cycle expansions is vanishing. For a > 0, Aq = /'(O) 1, 
the fixed point is marginally stable and is also customar- 
ily omitted from cycle expansions [^. The intermittent 
behavior arises from cycles containing long strings of O's 
which come close to the marginally stable fixed point. 



III. CYCLE EXPANSIONS 

Cycle expansion approaches to deterministic diffusion 
in one-dimensional maps were introduced in refs. [^ 18 
and in the Lorentz gas in ref. [ p^ . The dynamical zeta 
function formula ||] for the diffusion coefficient is 



D 



2 v^/ -I \fc "pi+"P2+---+"pfc 

^) IA„, A„„...A„, 



(4) 



where the sum is over nonempty distinct nonrepeating 
combinations of prime cycles, n is the lattice translation 
of a cycle, n is the period of the cycle, and A its stability. 
As the flow is conserved, the leading eigenvalue of the 
Frobenius-Perron operator equals unity, and the inverse 
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of the corresponding dynamical zeta function g| must 
vanish js) for z = 1: 



1/C(1) = 1 + 5](-l)V|A,,Ap, . . . ApJ = 



(5) 



For example, with the { — ,0,+} symbolic dynamics the 
cycle expansion up to topological length n = 2 equals 



1 



1 



A+A_ 



where we have omitted the cycle. 

Ideally the H — cycle is shadowed by the + and — cy- 
cles, so the last two terms are expected to approximately 
cancel. The cancellation is exact when a = and both 
terms equal 9. However, for other values of a shadow- 
ing may not lead to any significant cancellations. For 
example, for a — 1, A+A_ = 25 but A^ « 12. 

For the case a = the map f{x) is piecewise linear, 
A = 3" for all cycles, and in this case the stability and 
length ordering are equivalent. If the cycle is included, 
all n > 1 terms in (||) cancel, leading to D = 1/3. 

For a > 0, the dynamics is intermittent, and Ap does 
not necessarily grow exponentially with cycle length Hp. 
Furthermore, for a > the most stable orbits which con- 
tain long strings of O's are not shadowed by combinations 
of shorter cycles, as the cycle is not included. In fact, 
we can explicitly deduce the behavior of the most sta- 
ble cycles, those of the form O"-!- (an example is given in 
Fig. 1^) . These begin at initial point which we shall denote 
by Xn, slightly greater than zero. Many iterations of the 
function increase monotonically the value of /''{x) until 
it finally crosses over to the right branch and returns to 
the starting point. Inserting an extra in the symbolic 
dynamics has the effect of slightly decreasing the starting 
point, but the other cycle values are virtually unchanged, 
the new starting point is very close to the old one, and 
f{^n+i) = Xn is a good approximation for moderately 
large n. Thus for a > 

1 + 2(2a;„+i)" 



Xn+l 



to the leading order in x„ - This difference equation may 
be approximately solved as a power law, a;„ = •jn^^ giv- 
ing aS = 1, 2(27)" = S, or 

1 



" 2(2an)i/" ' 
The stability can be estimated in a similar fashion: 

^" nx„) = l + 2ia + l)i2xnr = 1+''^^ 



A„_i an 

Again putting A„ — pn'^ we obtain e = 1 + 1/a, 

A„ ^ ni+V" , 

confirmed by our numerical results. This power-law 
growth of A„ is in contrast to hyperbolic systems for 



which all cycles have stabilities which grow exponentially 
with the cycle length. 

Because these are the most stable cycles, they domi- 
nate cycle expansions at given n. Combinations of cycles 
which do not include a cycle with a string of almost n O's 
are highly suppressed. For example, two cycles with n/2 
O's and one other symbol each have a combined stability 



A2 



(nV4)'+'/" » A, 



for large n, again in contrast with hyperbolic systems, for 
which such shadowing combinations are of comparable 
magnitude. Hence we can estimate the convergence prop- 
erties of such cycle expansions by approximating them 
with the dominant 0"-f cycle family In this approxi- 
mation the flow conservation condition (ph 



i/c(i)-E' 



1 



,1+1/a 



(6) 



is approximately the Riemann zeta function C(l + 1/a), 
convergent for all a > 0, which is just as well. The de- 
nominator in the diffusion formula (^) appears whenever 
we calculate the time average of some quantity, and plays 
the role of a mean cycle period 



1 



(7) 



As a increases, the system spends more and more of 
its time near the marginally stable fixed point, and for 
a > 1, the system spends on average all of its time within 
an arbitrarily small neighborhood of this point, leading 
to a divergent mean cycle period (|^). The numerator of 
the expression for diffusion looks like the flow conserva- 
tion sum (||), with extra factors of rip. This factor is a 
number of order unity for the least unstable cycles, so 
the series converges. Thus the average which defines 
the diffusion coefficient undergoes a phase transition [T^ ] 
and equals zero for a > 1. This behavior is described as 
"weak" (0 < a < 1) or "strong" (a > 1) intermittency. 
Other averages may converge for different ranges of a. 

For a > 1 the diffusion is anomalous, with (n^) in- 
creasing more slowly than n. In the case at hand 




a < 1 
a = l 
a> 1 



and the leading behavior of a dynamical zeta function as 
z — > 1 yields |lq ] the exponent characterizing the sublin- 
ear diffusion 



n a < 1 

n/ Inn a — 1 
n^/" a > 1 



We are now in position to estimate and compare the 
rates of convergence of the topological and stability trun- 
cation approaches. As we have seen, the cycle expan- 
sions are dominated by terms of the form l/n"' where 
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7 = 1 + 1/a for the flow conservation condition, and 
7 = 1/a for the diffusion constant. Thus the error made 
by the topological truncation length after n terms is 
~ n}~'^ . In contrast, truncating by stability corresponds 
to an error of order of A^^/'^"+^) for the flow conserva- 
tion, and A^"^^^/("+^^ for the diffusion constant. For 
a — > 1 close to where the expansion diverges it may be 
advisable to improve the estimates by convergence accel- 
eration techniques. 

Having deflned the cycle expansions and analyzed their 
behavior in the intermittent case, we now pose the more 
pragmatic question: What is the optimal ordering in 
practice? In the approximation we have been using, with 
only a single family of cycles contributing, the ordering 
is self-evident. However, the full expansion is only con- 
ditionally convergent, and as we have no proof that the 
stability ordering yields the correct results, our justiflca- 
tion will come from heuristic arguments, together with 
the numerical results. 

The topological length cutoff corresponds to a com- 
plete partitioning of the phase space into 3" periodic 
point neighborhoods, irrespective of the relative sizes of 
these neighborhoods. What is the meaning of the sta- 
bility cutoff? A cycle expansion can be interpreted ||] 
as a partition of the dynamical phase space into neigh- 
borhoods of periodic points i £ p, each of size ^ 
A fixed stability cutoff A selects a uniform partition of 
the phase space into ~ A regions, each of period approx- 
imately N = In A, the time needed for a neighborhood 
of a hyperbolic orbit to spread across the entire system. 
Each prime cycle has about N periodic points, so the 
number of prime cycles up to a given stability grows as 
A/ In A. This estimate is known as the dynamics version 
of the prime number theorem for Axiom A systems, given 
in rcf . . We find that the estimate is valid numerically 
for nonhyperbolic systems as well, in the case at hand for 
all values of a > —1. 

The dramatic difference between the two approaches 
is the number of cycles required in each case. The num- 
ber of prime cycles up to a given length increases expo- 
nentially with the length, in our case as ~ S^/n. The 
number of prime cycles up to a given stability is grows as 
A/ In A. Superficially, the topological ordering requires 
an exponential number of cycles, while the stability or- 
dering requires only a power law. The issue is how small 
is the error for a given truncation. In the case of nice 
hyperbolic flows this error is superexponentially small, 
but for intermittent systems, the size of the error is not 
known, and we have to resort to numerics to estimate it. 



IV. NUMERICAL RESULTS 

For a simple one-dimensional map with a complete 
symbolic dynamics, such as the map (|^) studied here, 
almost any reasonable cycle finding method jsj should 
yield thousands of prime cycles. The accuracy of 1/A's 



that we calculate approaches the machine precision. We 
implement the stability ordering by noting that for this 
map any cycle containing an extra symbol is less stable 
than the preceding one. We recursively increment cycle 
lengths, starting with + and — , and stopping when the 
stability cutoff is reached. The stability ordering is fast, 
as cycles containing large numbers of O's which dominate 
the expansion appear only a few times in this enumera- 
tion. The distribution of cycles as a function of Up and 
Ap is shown in Fig. 0. 



12 
10 
8 

In A 6 

4 

2 

12 
10 
8 

In A 6 



I I 
■ I 



I i • 



a = -0.5 



2 4 6 8 10 12 14 16 18 20 



a = 




8 10 12 14 16 18 20 



FIG. 2. The distribution of cycle eigenvalues as a function 
of the topological cycle length: For a = 1, the intermittent 
case, a fixed topological length n cutofi' misses many of the 
least unstable (but longer period) cycles. 

In Fig. ^ we present a comparison of how well differ- 
ent truncations respect the flow conservation rule (|^). 
There is a definite improvement as the stability cutoff is 
increased, demonstrating convergence. For smaller val- 
ues of a there is a significant amount of scatter due to a 
small number of unbalanced shadowing terms which vary 
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rapidly with a, however the error is consistently small. 
For example, for a < 1 the error curve corresponding to 
A ^ 10^ cutoff always hes below e'^ « 0.0025 . 



ll,|C-'(0.1)| -6 




FIG. 3. Numerical test of the flow conservation for dif- 
ferent truncations. The logarithm of the magnitude of the 
zeta function (pH is plotted as a function of a for stability 
cutoffs A = 10 , 10'^, 10'', and 10^ (plusses, open triangles, 
squares and circles, respectively), as well as for the topologi- 
cal length cutoff A*' = 10 (filled circles). The solid curves for 
a > 0.5 are the behavior expected from the previous section. 
In (l.5A~"~^), for the cutoff values of A, with the constant 
1.5 chosen to fit the data. 

The filled circles in Fig. ^ are obtained by using all 
cycles with N < 10, corresponding to roughly the same 
computational effort as the stability cutoff of 10^. The 
error is smooth but large, comparable to the stability 
cutoff only at a = 0, the solvable case. 




FIG. 4. Diffusion constant (|3|) as function of a computed 
by 3 different methods. The dotted line is obtained by di- 
rect numerical evaluation of (m) for 3 x 10* random initial 
points, each evolved by 3 x 10 iterations of the map. The 
dot-dash line is the cycle expansion average truncated at 
topological length N — 10, corresponding to 9381 cycles with 
omitted. The expansion is accurate only near a — 0, where 
the shadowing is exact. The solid lines are cycle expansions 
(^) truncated according to stability, with cutoff A — 10^ for 
the uppermost line at a = 2, followed by 10^, 10* and 10^ for 
the other curves. The dashed line is an extrapolation of 
these curves, a — 1 corresponds to a phase transition point 
beyond which D = 0; all numerical methods fail here. 



The diffusion constant, evaluated using three different 
methods, is plotted as a function of a in Fig. ^. Each of 
the three methods used about half an hour of computing 
time for each value of a. There is no analytic expression 
for D in general, except for D = 1/3 at a = and D = 
for a > 1, as explained above. 

The cycle expansions (||) truncated at topological 
length A*" = 10 (dot-dash line) give very poor estimates 
of 13 as a increases away from zero into the intermit- 
tent regime, as many of the least unstable cycles are not 
included (Fig. ||). More surprisingly, even for a < 0, 
the hyperbolic phase dominated by exponentially many 
unstable cycles, the estimate is poor, this time in the op- 
posite direction, because there are many omitted cycles 
near the -I- and — cycles with large displacements and 
moderate stability, A ^ A". The stability of the + and 
— cycles is A-|_ = 3 + 2a which is relatively small in this 
case. 

The cycle expansions (jj) truncated according to sta- 
bility (solid lines) approach the direct simulation (dotted 
line) as the cutoff is increased from 10^ (upper curve 
at a = 2) to 10^. Since the number of prime cycles less 
than A is asymptotically equal to A/ In A, we expect that 
about 10^/ In 10^ « 8700 cycles would be needed, not far 
from the number of cycles actually found, which lies in 
the range 8871-10066 for all the a shown. 

Due to the phase transition all three methods are only 
logarithmically convergent to zero at a = 1. Direct sim- 
ulation cannot yield accurate estimates of D near this 
point; the stability and topological length truncations 
could easily be improved by using the analytic struc- 
ture of the 0"-f series j|] for the map at hand, but as 
we are unlikely to have this information available in a 
generic case, such improvements are outside the scope of 
the present investigation. 

As discussed in the previous section, we expect that 
the error should scale as a power of the cutoff stabil- 
ity, related to a, or exponentially with n for a cutoff of 
10". If the convergence of a sequence is exponential and 
reasonably smooth, it should be possible to extrapolate 
to the limit using a variety of convergence acceleration 
schemes. We have tried Aitken's 5^-process If 
t„ and tn+i are three consecutive terms of a sequence, an 
improved estimate is 



tr. 



{tn+1 — tnY 

+1 ~ 2t„ + t„_i 



(8) 



Using this formula for the cycle expansions computed to 
a stability of 10'^, 10'' and 10^ we obtain the dashed line 
in Fig. ^, which is comparable to the direct simulation, 
and much closer to the true value of D than the value 
obtained by the length cutoff. 

It should be noted that the estimate (||) only works if 
the sequence is relatively smooth. Cycle expansions with 
stability cutoff are not particularly smooth as a function 
of the cutoff, for example see Fig. 5 of ref. [|j . This is 
because at each stage a small number of shadowing com- 
binations are unbalanced by the cutoff, and the number 
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of such mismatches varies rapidly with the cutoff. In 
this case, the convergence over the range lO'^-lO^ is suffi- 
ciently smooth to use (||) , however a smaller spacing such 
as 2.5 X lO**, 5 x 10^, lO'^ is dominated by fluctuations. 



V. CONCLUSION 

Our results indicate that cycle expansions may be used 
to calculate averages for intermittent systems with accu- 
racy comparable to direct simulations, as long as a stabil- 
ity cutoff is used. Stability ordering has a great simplicity 
in that it requires no knowledge of the dynamics, except 
what is contained in a finite cycle set. We conclude with 
a few possibilities for future directions. 

First, it is good to see rapidly converging expansions, 
but another thing to have rigorous limits to guarantee 
convergence. Chaotic systems often behave more nicely 
than it is possible to prove, however it would be advan- 
tageous to extend the proofs of superexponential conver- 
gence for length ordered cycle expansions of analytic hy- 
perbolic systems to the stability ordered case, hopefully 
allowing a wider class of dynamical systems. 

The stability ordering exhibits imperfect shadowing, 
which can lead to scatter in the results, as observed in 
Fig. ^. One possible remedy to this problem might be 
to replace the factor l/H^ by /(H^)/]!^ where the 
smoothing function / moves continuously from 1 to 0. 
This must certainly improve the shadowing, but it re- 
quires that cycles be found up to the largest stability at 
which / is non-zero, without utilizing these cycles fully. 
For our diffusion coefficient calculations, the results are 
smooth enough to use Aitken's method, so additional 
smoothing is probably unnecessary. 

We also note, that because the number of cycles less 
than a given A is roughly A/ In A, independent of the 
dimension of the space, stability ordering should be ap- 
plicable to high dimensional systems. In particular, the 
detailed structure of the dynamics need not be known, 
only an algorithm for finding the cycles in the first place, 
for example tracing out a long trajectory and looking 
for near repeats, which are then refined by some form of 
Newton's method. 

Finally, it is not clear to what extent the stability cutoff 
approach is applicable to quantum systems. It is not as 
easy to estimate the rate of convergence in this case, even 
for nice hyperbolic ffows, because the terms are complex, 
and the alternation of the Maslov phase within families of 
cycles analogous to the 0"-|- family studied above crucial 
for quantum convergence can lead to large errors in the 
stability cutoff cycle expansion truncations 62] . 
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